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Abstract 

The critical thermodynamics of an MA^-component field model with cubic 
anisotropy relevant to the phase transitions in certain crystals with complicated 
ordering is studied within the four-loop e expansion using the minimal subtrac- 
tion scheme. Investigation of the global structure of RG flows for the physically 
significant cases M = 2, N = 2 and M = 2, N = 3 shows that the model has an 
anisotropic stable flxed point with new critical exponents. The critical dimen- 
sionality of the order parameter is proved to be equal to = 1.445(20), that 
is exactly half its counterpart in the real hypercubic model. 
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We study the critical behaviour of an MA^-component field model with cubic 
anisotropy having a number of interesting applications to phase transitions in three- 
dimensional simple and complicated systems. The effective Ginzburg-Landau Hamil- 
tonian of the model reads: 



-I N N N 
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where each vector field has M real components^ iff, i = 1, . . . , M and d = 4 — e 
is the spatial dimensionality. Here tjIq ~ (T — T^) and niQ, uq, vq are the "bare" mass 
and coupling constants, respectively. 

For M = N = 2 the Hamiltonian (|I]) describes the structural phase transition in 
Nb02 crystal and the antiferro magnetic phase transitions in TbAu2 and DyC2. Another 
physically important case M = 2, = 3 is relevant to the antiferromagnetic phase 
transitions in K2lrCl6, TbD2 and Nd materials The magnetic and structural phase 
transitions in a cubic crystal are governed by the model (|l]) at M = 1 and = 3 
0. In the replica limit — (M = 1) the Hamiltonian (|I|) is known to determine 
the critical properties of weakly disordered quenched systems undergoing second-order 
phase transitions with a specific set of critical exponents [Q. Finally, the case M = 1 
and N ^ oo corresponds to the Ising model with equilibrium magnetic impurities . 
In this limit the Ising critical exponents take Fisher- renormalized values |0. Since 
static critical phenomena in a cubic crystal as well as in randomly diluted Ising spin 
systems are well understood [0, ||, ^ |10], we will focus here on the critical behaviour 
of the above mentioned multisublattice antiferromagnets. This is the case of M = 2 
and A^ = 2, A^ = 3 in the fiuctuation Hamiltonian (|l|). 

For the first time, magnetic and structural phase transitions in crystals with com- 
plicated ordering described by the model (U) were studied by Mukamel and Krinsky 
within the lowest orders in e: [0. A stable fixed point ("unique" point), different from 
the isotropic (O(MA^)-symmetric) or the Bose (O(MA^)-symmetric, MN = 2) one[|, 
was predicted in d = 3. The point was shown to determine a new universality class 
with a specific set of critical exponents. However, for the physically important case 
M = N = 2, the critical exponents of the unique fixed point turned out to be the same 
as those of the 0(4)-symmetric one within the two- loop approximation. 

Later an alternative analysis of critical behaviour of the model, the RG approach 
in fixed dimension, was carried out within the two- and three-loop approximations 



1^ . Those investigations gave the same qualitative predictions: the unique stable 
fixed point does exist on the three-dimensional RG fiow diagram. However, the critical 
exponents computed at this point with the use of different resummation procedures 



-'^For M = 1 and M — 2 the model (0) is merely the iV-component cubic model determined either 

by the real or by the complex order parameter field, respectively. 

^For a six-loop study of the critical behaviour of the random Ising model see [|2| 

■^Note that the isotropic or the Heisenberg fixed point can be determined as the point having the 

coordinates u > 0, u = on the RG flow diagram, whereas for the Bose fixed point one has u — 0, 

v>0. 



1 



proved to be close to those of the Bose fixed point rather than the isotropic one. It 
was also shown that the unique and the Bose fixed points are very close to each other 
on the diagram of RG fiows, so that they may interchange their stability in the next 



orders of RG approximations [Tl| 



Recently, the critical properties of the model were analyzed in third order in e 
p!5| , [1^. Investigation of the fixed point stability and calculation of the critical di- 
mensionality Nc of the order parameter, separating two different regimes of critical 
behaviour, confirmed that the model (|l]) possesses the anisotropic (complex cubic; 
M = 2 and -u 7^ 0, 7^ 0) stable fixed point for N = 2 and 3. The realistic critical 



exponent estimates for the unique stable fixed point were obtained in |T6[ using the 



summation technique of [T^ . The values appeared to be close to those of the isotropic 
point in contradiction to the numerical results given by RG directly in three dimensions 



p3| , Such a distinction may be accounted for by the too short three- loop e series 
used. 

It is worthy of note that the existence of an anisotropic stable fixed point on the 
three-dimensional RG fiow diagram contradicts the nonperturbative considerations 



T8|| . Indeed, according to those considerations the only stable fixed point in three 
dimensions may be the Bose one and it is that point which governs the critical thermo- 
dynamics in the phase transitions of interest. The point is that the model (|l]) describes 
N interacting Bose systems. As was shown by Sak , the interaction term can be rep- 
resented as the product of the energy operators of various two-component subsystems. 
It was also found that one (the smallest) of the eigenvalue exponents characterizing 
the evolution of this term under the renormalization group in a neighborhood of the 
Bose fixed point is proportional to the specific heat exponent a. Since a is believed to 
be negative at that point, and that is confirmed experimentally pO[ and theoretically 



21[| , the interaction is irrelevant. Consequently, the Bose fixed point should be stable 
in three dimensions. However, the RG approach, applied to the model (|I|), has not 
yet confirmed this conclusion. It is therefore highly desirable to extend already known 
e expansions for the stability matrix eigenvalues, critical exponents and the critical 
dimensionality in order to apply more sophisticated resummation techniques to longer 
expansions. 

Therefore, the aim of this Letter is to extend the existing three-loop e expansions 
of the model to the four-loop order and to study more carefully the predictions of the 
RG method regarding the critical phenomena in the substances of interest. Namely, on 
the basis of the four-loop expansions for the RG functions obtained using dimensional 



regularization and the minimal subtraction scheme p2| , p^] , we analyze the stability 
of the fixed points and calculate the critical dimensionality of the order parameter. 
We show that the anisotropic stable fixed point does exist on the thre-dimensional RG 
fiow diagram. For this point, we give more accurate critical exponent estimates, in 
comparison with the previous three-loop results by applying a new summation 
approach [17| to the four- loop series. 
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The four-loop e expansions for the /^-functions of the model are as follows 



eu — u 
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X (3528Ar + 21240) - 480C(5) (lOAT + 287) + 61 A^^ - 5173Ar - 66764) - 

X (1800C(3)(A^ + 62) - 144C(4)(8A^ + 203) + 172800C(5) + 56A^ + 93701) 

- 4mv^(5090C(3) - 1296C(4) + 7600C(5) + 4503) + ^(-8224C(3) + 1920C(4) 

- 12160C(5) - 7975) 

where C(3),C(4), and C(5) are the Riemann C functions. 

Prom the system of equations /?„(«*, f*) = and /3„(w*,f*) = 0, we found formal 
series for the four fixed points: trivial Gaussian and nontrivial isotropic, Bose, and 
complex cubic. Then we calculated the eigenvalues of the stability matrix 

C d(3uiu,v) d(3uiu,v) \ 
du , dp ^ \ 

df3^iu,v) df5^(u,v) 
du dv ' 

taken at the most intriguing Bose and complex cubic fixed points. The corresponding 



numerical estimates are obtained using an approach based on the Borel transformation 



Fie;a,b) = ±MX) T e"^ (-)'<-) t^^^^ (4) 



modified with a conformal mapping z = "^^^^j^^ P^ , which does not require the knowl- 



edge of the exact asymptotic high-order behaviour of the series 17|. The parameter 



A is chosen from the condition of the most rapid convergence of series (^), that is by 
minimizing the quantity |1 — p^^^^^^^l, where / is the step of truncation and Fi{e; a, b) 
is the /-partial sum for F{e; a,b). If the real parts of both eigenvalues are negative, 
the associated fixed point is infrared stable and the critical behaviour of experimental 
systems undergoing second-order transitions is determined only by that stable point. 
For the Bose and the complex cubic fixed points our numerical results are presented 
in table 1. It is seen that the complex cubic fixed point is absolutely stable in d = 3 
{e = 1), while the Bose point appears to be of the "saddle" typeQ. However, the uj2 of 
either point is very small at the four-loop level, thus implying that these points may 
swap their stability in the next order of the RG approximation. 

In addition to the eigenvalues, we calculated the critical dimensionality A*"^ of the 
order parameter from the condition of vanishing (jj2 for the complex cubic fixed point. 
The four-loop expansion is 



A^f = 2 - £ + -- 6C(3) -I e' + -^ 45C(3) + 135C(4) - 600C(5) - 1 e\ (5) 
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24 



1 



144 
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Instead of processing this expression numerically, we established the exact relation 
= |iV^, which is independent on the order of approximation used^. In fact, the 
critical dimensionality for the complex cubic model is determined as that value of 
A^*^ at which the complex cubic fixed point coincides with the isotropic one. The same 
assertion holds for the cubic model with the real A^^-component order parameter. So, 
because of the relation 0(2A^'^) = 0{N^), the relation 2A^f = A^^^ should hold too. 

The five-loop e-expansion for was recently obtained in Ref. 0. Resummation 
of that series gave the estimate A"^ = 2.894(40) 0. Therefore we conclude that 
= 1.447(20) from the five- loops. Practically the same estimate = 1.435(25) 
follows from a constrained analysis of A^^ taking into account A"^ = 2 in two dimensions 



T0|| . From the recent pseudo-5 expansion analysis of the real hypercubic model [|ri|] 
one can extract A^^ = 1.431(3). However the most accurate estimate A"^ = 1.445(20) 
results from the value = 2.89(4) obtained on the basis of the numerical analysis 
of the four-loop [^] and the six-loop three-dimensional RG expansions for the f3- 
functions of the real hypercubic model. Since < 2, the phase transitions in the 
Nb02 crystal and in the antiferromagnets TbAu2, DyC2, K2lrCl6, TbD2 and Nd are of 



^We say the fixed point is of the "saddle" type provided their eigenvalue exponents uJi and llJ2 are 
of opposite sings in the {u, v) plane. 

^Here and are the critical (marginal) spin dimensionalities in the complex (A/ — 2) and 
in the real (M = 1) hypercubic model, respectively. 
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second order and their critical thermodynamics should be controlled by the complex 
cubic fixed point with a specific set of critical exponents, in the frame of the given 
approximation. The corresponding four-loop critical exponent estimates are displayed 
in table 2. The critical exponent estimates obtained for the isotropic and the Bose 
fixed points are also presented in the table, for comparison. 

In conclusion, the four-loop e-expansion analysis of the generalized MA^-component 
Ginzburg-Landau model with cubic anisotropy describing phase transitions in certain 
real antiferromagnets with complicated ordering has been carried out with the use of 
the minimal subtraction scheme. Investigation of the global structure of RG fiows for 
the physically significant cases M = 2, = 2 and M = 2, = 3 has shown that 
the complex cubic rather than the Bose fixed point is absolutely stable in three dimen- 
sions. The critical dimensionality A*"^ = 1.445(20) of the order parameter obtained 
from six loops has confirmed this conclusion. For the stable complex cubic fixed point, 
reasonable estimates of critical exponents were obtained using the Borel summation 
technique in combination with a conformal mapping. For the structural phase tran- 
sition in Nb02 and for the antiferromagnetic phase transitions in TbAu2 and DyC2, 
they were shown to be close to the critical exponents of the 0(4)-symmetric model. 
In contrast to this, the critical exponents for the antiferromagnetic phase transitions 
in K2lrCl6, TbD2 and Nd turned out to be close to the Bose ones. Although our re- 
sults seem to be substantially consistent with other predictions of the RG approach, 
still there is a definite contradiction with the general nonperturbative theoretical argu- 
ments []TB| mentioned in the introduction. One can hope, however, that the five-loop 



contributions being taken into account will eliminate this controversy. Indeed, the 
present calculations have shown that, although the complex cubic fixed point, rather 
than the Bose one, is stable at the four-loop level, the eigenvalues 002 of both fixed 
points are very small. Therefore the situation is very close to marginal, and the points 
might change their stability to the opposite in the next order of perturbation theory, 
so the Bose point would turn out to be stable. There is a hope that comparison of 
the critical exponent values obtained theoretically for different fixed points with those 
values determined from experiments or, probably, from Monte Carlo simulations would 
indicate which fixed point is really stable in three dimensions. Finally, it would also be 
desirable to investigate certain universal amplitude ratios of the model because they 
vary much more among different universality classes than exponents do and might be 
more effective as a diagnostic tool. 

We are grateful to Professor M. Henkel for helpful remarks and to Dr. E. Bla- 



goeva for communicating her results mentioned in [15|. One of the authors (K.B.V.) 
acknowledges useful discussions with Dr. B. N. Shalaev. This work was supported by 
the Russian Foundation for Basic Research, grant no 01-02-17048, and by the Ministry 
of Education of Russian Federation, grant no EOO-3. 2-132. 

Note added. More recently, the six-loop RG functions of the model of interest have been 



calculated within the alternative RG approach in three dimensions |T^, 25 1. Although 



the authors argue the global stability of the Bose fixed point, the numerical estimate of 
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Table I: Eigenvalue exponents estimates obtained for the Bose (BFP) and the complex 
cubic (CCFP) fixed points at = 2 and = 3 within the four-loop approximation in 
e (e = 1) using Borel transformation with a conformal mapping. 



Type of 
fixed point 


A^ = 2 


A^ = 3 




Ll)2 




Ll)2 


BFP 


-0.395(25) 


0.004(5) 


-0.395(25) 


0.004(5) 


CCFP 


-0.392(30) 


-0.029(20) 


-0.400(30) 


-0.015(6) 



Table II: Critical exponents calculated for the isotropic (IFP), the Bose (BFP), and 
the complex cubic (CCFP) fixed points at A^ = 2 and A^ = 3 within the four-loop 
approximation in e (e = 1) using Borel transformation with a conformal mapping. 



Type of 
fixed point 


A^ = 2 


A^ = 3 


7] 


V 


7 


ri 


V 


7 


IFP 


0.0343(20) 


0.725(15) 


1.429(20) 


0.0317(10) 


0.775(15) 


1.524(25) 


BFP 


0.0348(10) 


0.664(7) 


1.309(10) 


0.0348(10) 


0.664(7) 


1.309(10) 


CCFP 


0.0343(20) 


0.715(10) 


1.404(25) 


0.0345(15) 


0.702(10) 


1.390(25) 



the smallest stability matrix eigenvalue of the Bose fixed point, which governs the RG 
fiows near this point, appears to be very small, {jJ2 = —0.007(8) [^, and the apparent 
accuracy of the analysis does not exclude the opposite sing for UJ2. This result agrees 
well with our conclusions. 
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